Interplay between Plasmon Luminescence and Vibrationally Resolved Molecular 
Luminescence Induced by Scanning Tunneling Microscopy 
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Effects of coupling between a molecular exciton and a surface plasmon (exciton-plasmon coupling) 
on the luminescence properties of the molecule and the surface plasmons are investigated using the 
nonequilibrium Green's function method. Molecular absorption and enhancement by molecular elec- 
tronic and vibrational modes (molecular modes) lead to dip and peak structures in the luminescence 
spectra of the surface plasmons. These structures will correspond to the peak and dip structures 
observed in a recent experiment. We found that in addition to the molecular dynamics, the re- 
absorption by the surface plasmons plays important roles in determining the luminescence spectral 
profiles. 

PACS numbers: 73.20.Hb, 73.20. Mf, 33.50.Dq 



Fundamental understanding of the role of surface plas- 
mons in the molecular luminescence process is currently 
a topic of great interest owing to its wide applicability 
in nanoscale spectroscopy, biosensors, light-emitting de- 
vices, etc [1|-|6[. One of the most promising approaches 
for deepening our knowledge is to investigate light emis- 
sion induced by tunneling current of a scanning tunneling 
microscope (STM). STM-induced light emission (STM- 
LE) spectroscopy has a unique advantage of having the 
potential to unveil optical properties with submolecular 
spatial resolution. Moreover, it can remove any complex- 
ity resulting from disturbance of the incident light. 

STM-LE from a metal surface was first reported by 
the group of Gimzewski in 1988 (7|. Now, it is well con- 
firmed that light emission from clean metal surfaces is 
due to the radiative decay of surface plasmons localized 
near the tip-sample gap region When a molecule 

is inserted into the gap region, there are two radiative 
processes, i.e., luminescence from the surface plasmons 
and the molecule. If the molecule is directly adsorbed on 
a metal substrate, charge and energy transfer between 
them lead to the quenching of molecular luminescence. 
Thus, the molecule acts as a spacer between the tip and 
the metal substrate, which slightly modifies the plasmon- 
mediated emission [l^]- If the molecule is electrically 



decoupled from the metal substrate by dielectric films or 
molecular multilayers, intrinsic molecular luminescence, 
which is associated with the electronic and vibrational 
transitions in the molecule, can be observed 13- 16|. 

Since close proximity of the tip to the metal substrate 
induces an intense electromagnetic field generated by the 
surface plasmons, effects of the interaction between the 
intense electromagnetic field and the transition moments 
for the molecular excitations and de-excitations are ex- 
pected to occur. Due to the fact that the molecular lu- 
minescence can be in the same energy window as the 
plasmon-mediated emission, light emission can be in- 
duced by the two kinds of radiative processes simultane- 
ously. Recent experiments have shown that the plasmon- 



mediated emission overlapping with the molecular ab- 
sorption enhances molecular luminescence intensity 
ITgi l . Selective enhancement of luminescence intensity by 
spectrally tuning the surface plasmon mode to match 
with a particular transition in the molecule has been 
observed [20| . Moreover, fascinating new phenomena, 
i.e., luminescence associated with radiative decay from 
higher vibrational levels for the excited electronic state 
of the molecule (hot luminescence) and luminescence at 
energies beyond the applied bias voltage (upconverted 
luminescence) have also been observed [20|. 

There are two dominant mechanisms involved in the 
molecular excitations in STM-LE from molecules on 
metal substrates. One is the excitation by the injection 
of electrons and/or holes from the electrodes. The other 
is the excitation by the absorption of surface plasmons 
that are excited by the inelastic tunneling between the tip 
and the substrate. Recent experiments suggest the im- 



portance of the plasmon-mediated excitations |17H21|. 
In previous theoretical studies, molecular dynamics have 
been investigated within the framework of the one-body 
problem, where the surface plasmon is treated as a clas- 
sical electric field fH-flll • Tian et al. [H, [ll] used the 
plasmon-mediated excitation mechanism to explain some 
spectral features observed in Ref. IIS- The results indi- 
cate that the molecular dynamics induced by the surface 
plasmons play essential roles. 

Direct experimental evidences of the plasmon- 
mediated excitations of the molecules have been obtained 
in luminescence spectra acquired with molecule-covered 
tips over clean metal surface |25|]. Although no electron 
tunneling to the molecules takes place, the observed spec- 
tra that can be considered as the luminescence spectra of 
the surface plasmons show peak and dip structures, some 
of which match peaks in molecular luminescence and ab- 
sorption spectra. Thus, the molecular dynamics includ- 
ing the molecular luminescence and absorption have an 
impact on luminescence-spectral profiles of the surface 
plasmons. To understand this from a microscopic point 



2 



plasmon 

■^f\l\r* 




Energy 



Q 




Qo 



FIG. 1. (Color online) Schematic energy diagram of exci- 
tation in the absorbed molecular layer whose excitation en- 
ergy is close to the surface plasmon mode. Horizontal lines 
in each parabola denote vibrational sublevels where \g) and 
|e) denote the ground and first-excited electronic states, re- 
spectively. The variable Q denotes the vibrational coordinate 
whose equilibrium position is displaced by Qo upon the cre- 
ation of the molecular exciton. 



of view, there is a need to investigate the interplay be- 
tween the dynamics of the molecule and the surface plas- 
mons within the framework of the quantum many-body 
theory. 

In this study, we investigate the effects of coupling 
between a molecular exciton and the surface plasmon 
(exciton-plasmon coupling) with the aid of the nonequi- 
librium Green's function method [26| as well as coupling 
between electronic and vibrational degrees of freedom 
on the molecule (exciton- vibron coupling). Our results 
show that the enhancement and suppression of the lu- 
minescence intensities of the molecule are due to the 
enhancement by the surface plasmons and due to the 
re-absorption of the surface plasmons by the molecule, 
respectively. Absorption and enhancement by molec- 
ular electronic and vibrational modes lead to dip and 
peak structures in the luminescence spectra of the sur- 
face plasmons. In addition to the molecular dynamics, it 
was found that the re-absorption by the surface plasmons 
leads to a dent structure in the luminescence spectra of 
the surface plasmons. The results also show evidence 
that vibrational excitations assist the occurrence of the 
upconverted luminescence. 

A schematic illustration of our model is shown in Fig. 
[TJ The electronic degrees of freedom of the molecule are 
modeled as a system with the ground electronic state \g) 
and the first-excited electronic state |e). The molecu- 
lar exciton interacts with the molecular vibrations and 
with the surface plasmons. For simplicity, the molecular 
vibrations are approximated by a single harmonic oscil- 
lator [ij] . The surface plasmon modes localized near the 
tip-sample gap re gion are approximated by a single en- 



ergy mode (20|, |2J, [27| . The creation (annihilation) of the 



tem is 

13 

+AIQbdU + V (ad^ + H.c.) + ^ UfsQtQp, (1) 

13 

where and d are creation and annihilation operators 
for the molecular exciton with energy eex- Operators 
6^ and b are creation and annihilation operators for a 
molecular vibrational mode with energy hojQ, and a 
are for a surface plasmon mode with energy huip, bjj and 
bp are for a phonon mode in the thermal phonon bath, 
Qb — b + b''' and Qp — bp + fejj. The energy parameters 
M, V and Up correspond to the exciton-vibron coupling, 
the exciton-plasmon coupling, and the coupling between 
the molecular vibrational mode and a phonon mode in 
the thermal phonon bath. 

By applying the canonical (Lang-Firsov) transforma- 
tion H becomes 



H =eexd^d + hojob^b + hwpo) a + hwpb^^pbp 

f3 

+V [ad^X^ + H.c.) + UpQbQp, 



(2) 



molecular exciton is induced by the absorption (emission) 
of the surface plasmon 25] . The Hamiltonian of the sys- 



where X = cxp [—A (b^ ^ b)] , e^x = £ex ~ M^/ (hujo) and 
A = M/{hujo). 

The luminescence spectra of the molecule and the sur- 
face plasmons are obtained by the relations, Bl{oj) = 
-3L<(w)/7r and Bp{uj) = -!3P<(w)/7r Here L< 

and P< are the lesser projections of the Green's function 
of the molecular exciton and the surface plasmon on the 
Keldysh contour (L and P), respectively. 

To investigate the excitation properties and analyze 
the structures in the luminescence spectra, the spec- 
tral functions of the molecule and the surface plasmons, 
which correspond to their absorption spectra, are cal- 
culated by the relations, Al{u!) = — 3L''(a;)/7r and 
Ap{uj) = —^P^ {uj)/tt, where and P'' are the retarded 
projections of L and P. 

To obtain the spectral functions and the luminescence 
spectra, it is needed to calculate L and P, which are 
defined as 

L(r,r') ^ l(Tc{d(r)X(r)dt(r')Xt(r')})^, (3) 

Pir,r')^^{Tc{a{r)aHr')})^, (4) 

where (• • • ) denotes statistical average in representation 
by system evolution for H. The variable r is the Keldysh 
contour time variable. The operator Tc denotes the time 
ordering along the Keldysh contour. 

To calculate the Green's function L containing the 
exciton-plasmon coupling we take into account a 
random-phase-approximation (RPA)-type diagrammatic 
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series. This approximation is expected to hold for M ^> 
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V''^ [30|. The integral equation for L is given by 
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L(T,T')=ib(T,T') 
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iei(r,r') = ^(TcMr)dt(r')})^, 
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X(r,r') = (Tc{X(r)Xt(r')})^, 


(8) 
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and is the plasmon Green's function for 


y = 0. 






Here, we employ an approximation of decoupling elec- 
tronic and vibrational parts in Lb, and calculate Ld and 
K. The integral equation for Lei is given by 

Lei(r,T')-L(")(r,r') 

dTldT2L(")(T,Tl)I](Tl,r2)Lel(r2,T'), (9) 
S(ti,T2) = |ypp(0)(Tl,T2)if(T2,Tl), (10) 

where L'"' is Ld for V = 0. The integral equation for P 
is given by 

P(t, t')= 

+ / dridT2P(")(T,Tl)|ypL(Tl,T2)F(r2,T').(ll) 



We express the correlation function K in terms of the 
vibrational Green's function L>, which is defined by 

D{r,r') = ^{Tc{P,iT)P,{r')})^ (12) 

where Pb — ~i (b — 6^) [l^. The self-energy for D can 
be expressed by using P and L. 

By assuming the condition of stationary current, the 
distribution function A'pi of the surface plasmons excited 
by inelastic tunneling between the tip and the substrate 
is given by 



iVpi(w) 




eVbias 



\f^\ < leVbiasI 

others 



(13) 



where Tpi is a coefficient related to the current ampli- 
tude due to the inelastic tunneling [3l|, e is the elemen- 
tary charge, and 14ias is the bias voltage. The lesser and 
greater projections of P(°^ are obtained through the rela- 
tions, P(0)'<(a;) = 2iAfpi(w)SP(°)''-(a;) and P(°)'>(a;) = 
2i [I + 7Vpi(w)] 9P(°)''', where P(°)''' is the retarded pro- 
jection of P^°^ [2^. The equations mentioned above are 
self-consistently solved. After convergence, the spectral 
functions and the luminescence spectra are calculated. 

The parameters given correspond to the experiment 
on STM-LE from tetraphenylporphyrin(TPP) molecules 
on the gold surface [20[: iex = 1-89 eV, hioo = 0.16 eV 
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FIG. 2. (Color online) Spectral functions of the molecular 
exciton Al for (a) V = eV, (b) V = 0.05 eV, (c) V = 
0.10 eV and (d) V = 0.20 eV, and the surface plasmon Ap 
for {e)V = eV, (f) V = 0.05 eV, (g) V = 0.10 eV and (h) 
V = 0.20 eV. Red solid and green dashed lines show functions 
for the bias voltage Vbias = 1.8 V and 2.5 V, respectively. 



and hojp = 2.05 eV. The statistical average is taken for 
temperature T = 80 K 0|. The square of A has been 
reported to be 0.61 on the basis of first-principles cal- 
241 . The parameter Up is given so that the 



culations 

molecular vibrational lifetime due to the coupling to the 
thermal phonon bath is 13 ps A Markovian decay 
is assumed for the surface plasmon so that the plasmon 
lifetime for y = eV becomes 4.7 fs [10, S!]. The co- 
efficient Tpi is set to 10^* to give the tunneling current 
It = 200 pA. 

Figure [2] shows the calculated spectral functions Al 
and Ap. The peak positions for both Al and Ap are 
shifted due to the exciton-plasmon coupling V , remark- 
ably for the lowest one upon excluding the small structure 
near 1.7 eV. It is confirmed that the amount of the shift 
is roughly scaled by in the range from V = 0.05 cV 
to V — 1 eV. These peaks include the contribution from 
the processes of electronic transitions accompanied by vi- 
brational state transitions. The amplitude of the asym- 
metric structure in the vicinity of 1.7 eV increases as the 
bias voltage Vbias increases. These structures originate 
from a process which involves energy transfer from the 
vibrational state. 

The luminescence spectra Bl and Bp are shown in 
Fig. [3l At Vbias = 1.8 V, the luminescence intensity of 
the molecule [Fig. ^a)] is lower than the surface plas- 
mons [Fig. [Slje)]. While at Vbias > 1-85 V, the lumines- 
cence intensity of the molecule becomes higher than the 
surface plasmons. In addition, at Vbias = 2.5 V, the lu- 
minescence spectra Bl [red sohd fine in Fig. ^d)] shows 
peak structures near 2.05 eV and 2.21 eV. The position 
of these peak structures correspond to the energies of 
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FIG. 3. (Color online) Luminescence spectra of the molecule 
Bl for the bias voltage (a) Vbias = 1.8 V, (b) Vbias = 1.85 V, 
(c) Vbias = 2.0 V and (d) Vbias = 2.5 V, and the surface 
plasmon Bp for the bias voltage (e) Vbias ~ 1.8 V, (f) Vbias ~ 
1.85 V, (g) Vbias = 2.0 V and (g) Kias = 2.5 V. Red solid 
and green dashed lines in Sl are luminescence spectra for L 
and Li,, respectively. Red solid and green dashed lines in Bp 
are luminescence spectra for P and P'^-*, respectively. The 
exciton-plasmon coupling is V = 0.10 eV. 



transitions from the higher vibrational levels for the ex- 
cited electronic state to the ground electronic state of the 
molecule. Therefore, these peak structures are attributed 
to the hot luminescence, which have been observed in the 
experiment [20| . The luminescence intensity of the peak 
near 1.89 eV is suppressed. In this energy range, the 
high intensity peak is found in the spectral function of 
the molecule [Fig. [2{c)]. Thus, the suppression of the lu- 
minescence intensity is attributed to the re-absorption of 
the surface plasmons by the molecule (re-absorption pro- 
cess). The re-absorption processes are classified into two 
kinds: one is accompanied by the vibrational excitations 
while the other is not (inelastic and elastic processes, re- 
spectively). From the spectral function of the molecule, 
it is confirmed that for < 1, the elastic process is 
dominant over the inelastic processes. 

In the luminescence spectra of the surface plasmons 
Bp, complicated peak and dip structures appear at 
Fbias > 1-85 V [red solid line in Figs. [Stg), [^h)]. The 
position of the peak structures in the energy range lower 
than 1.75 eV (near 1.70, 1.54 and 1.38 eV) correspond to 
the peaks in the luminescence spectra of the molecule Bl 
[red solid line in Figs. [3Uc),[Hd)]. These peak structures 
are, therefore, due to the enhancement by the molecu- 
lar electronic and vibrational modes (molecular modes). 
The position of the dip structures in the energy range 
higher than 1.8 eV in Bp (near 1.86, 2.02 and 2.18 eV) 
correspond to the peaks in the spectral function of the 
molecule A^. These dip structures are attributed to the 




1.85 1.90 1.95 
Energy (eV) 




1.85 1.90 1.95 
Energy (eV) 



FIG. 4. (Color online) Luminescence spectra of (a) the 
molecule Bl and (b) the surface plasmon Bp at the bias 
voltage Vbias ~ 1.8 V. Red solid and green dashed lines show 
spectra for vibrations in the nonequilibrium and thermal equi- 
librium states, respectively. The exciton-plasmon coupling is 
V = 0.10 eV. 



molecular absorption. These peak and dip structures 
will correspond to the peak and dip structures observed 
in the recent experiment psf . In addition, our results 
show that a dent structure appears near the energy of 
the surface plasmon mode (2.05 eV) in Bp [Fig. |3l[h)], 
although Al has its maximum intensities near 1.86 eV. 
It is found that the re-absorption by the surface plas- 
mons leads to the dent structure in the energy range near 
2.05 eV, where the spectral function of the surface plas- 
mons Ap has its maximum intensity. Thus, we insist 
that in addition to molecular luminescence and absorp- 
tion, the re-absorption by the surface plasmons play cru- 
cial roles in determining the luminescence spectral pro- 
files of the surface plasmons. We expect that these results 
can be verified by experiments, for example by a com- 
parison between the luminescence spectra acquired with 
a molecule-covered tip and with a clean metallic tip on 
a clean metal surface, which correspond to the red solid 
and green dashed lines in Fig. Elh), respectively 25|. 

The upconverted luminescence can be seen in Fig. 
m Its contribution can be confirmed in comparison 
with the vibrational state in thermal equilibrium [green 
dashed line in Fig. |3] , where the molecular vibration is 
distributed according to the Bose distribution function 
at T = 80 K and is therefore almost in the ground 
state. The vibrational excitations are accompanying 
the electronic transitions, which occur via the excitation 
channels resulting from the exciton-plasmon coupling V. 
The vibrational occupation number {h%)fj are scaled by 
\'^\V\'^nexM, where riex is the population of the molecu- 
lar exciton, which is given by riex — (d^d)^, and 7 rep- 
resents the vibrational damping for y = eV. Thus 
the vibrational excitations assist the occurrence of the 
upconverted luminescence. 

In conclusion, the luminescence properties of the 
molecule and the surface plasmons can be strongly in- 
fluenced by the interplay between their dynamics result- 
ing from the exciton-plasmon coupling. The lumines- 
cence spectra of the molecule are modified due to the 
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re- absorption of the surface plasmons by the molecule. 
The peak structures in Bp arise due to the enhancement 
by the molecular modes and the dip structures due to the 
molecular absorption. It is found that the re-absorption 
by the surface plasmons leads to the dent structure in 
Bp. Thus, in addition to the molecular dynamics, the 
re-absorption by the surface plasmons plays crucial roles 
in determining the spectral profiles of Bp. Moreover, our 
calculations have reproduced the hot luminescence and 
the upconverted luminescence assisted by the molecular 
vibrations. 
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